Merging KK-means with hierarchical clustering for identifying general-shaped groups

Clustering partitions a dataset such that observations placed together in a group are similar but different from those in other groups. Hierarchical and $K$-means clustering are two approaches but have different strengths and weaknesses. For instance, hierarchical clustering identifies groups in a tree-like structure but suffers from computational complexity in large datasets while $K$-means clustering is efficient but designed to identify homogeneous spherically-shaped clusters. We present a hybrid non-parametric clustering approach that amalgamates the two methods to identify general-shaped clusters and that can be applied to larger datasets. Specifically, we first partition the dataset into spherical groups using $K$-means. We next merge these groups using hierarchical methods with a data-driven distance measure as a stopping criterion. Our proposal has the potential to reveal groups with general shapes and structure in a dataset. We demonstrate good performance on several simulated and real datasets.Comment: 16 pages, 1 table, 9 figures; accepted for publication in Sta

Genome-Wide Association Study and Subsequent Exclusion of ATCAY as a Candidate Gene Involved in Equine Neuroaxonal Dystrophy Using Two Animal Models.

Equine neuroaxonal dystrophy/equine degenerative myeloencephalopathy (eNAD/EDM) is an inherited neurodegenerative disorder of unknown etiology. Clinical signs of neurological deficits develop within the first year of life in vitamin E (vitE) deficient horses. A genome-wide association study (GWAS) was carried out using 670,000 SNP markers in 27 case and 42 control Quarter Horses. Two markers, encompassing a 2.5 Mb region on ECA7, were associated with the phenotype (p = 2.05 × 10-7 and 4.72 × 10-6). Within this region, caytaxin (ATCAY) was identified as a candidate gene due to its known role in Cayman Ataxia and ataxic/dystonic phenotypes in mouse models. Whole-genome sequence data in four eNAD/EDM and five unaffected horses identified 199 associated variants within the ECA7 region. MassARRAY® genotyping was performed on these variants within the GWAS population. The three variants within ATCAY were not concordant with the disease phenotype. No difference in expression or alternative splicing was identified using qRT-PCR in brainstem across the ATCAY transcript. Atcayji-hes mice were then used to conduct functional analysis in a second animal model. Histologic lesions were not identified in the central nervous system of Atcayji-hes mice. Additionally, supplementation of homozygous Atcayji-hes mice with 600 IU/day of dl-α-tocopheryl acetate (vitE) during gestation, lactation, and adulthood did not improve the phenotype. ATCAY has therefore been excluded as a candidate gene for eNAD/EDM

New species longevity record for the northern quahog (=hard clam), Mercenaria mercenaria

Author Posting. © National Shellfisheries Association, 2011. This article is posted here by permission of National Shellfisheries Association for personal use, not for redistribution. The definitive version was published in Journal of Shellfish Research 30 (2011): 35-38, doi:10.2983/035.030.0106.Twenty-two large shells (>90 mm shell height) from a sample of live collected hard shell clams, Mercenaria mercenaria, from Buzzards Bay, Woods Hole, Cape Cod, MA, were subjected to sclerochronological analysis. Annually resolved growth lines in the hinge region and margin of the shell were identified and counted; the age of the oldest clam shell was determined to be at least 106 y. This age represents a considerable increase in the known maximum life span for M. mercenaria, more than doubling the maximum recorded life span of the species (46 y). More than 85% of the clam shells aged had more than 46 annual increments, the previous known maximum life span for the species. In this article we present growth rate and growth performance indicators (the overall growth performance and phi prime) for this record-breaking population of M. mercenaria. Recently discovered models of aging require accurate age records and growth parameters for bivalve populations if they are to be utilized to their full potential.This work was supported by grants from the American Diabetes Association (to Z. U.), American Federation for Aging Research (to A. C.), the University of Oklahoma College of Medicine Alumni Association (to A. C.), the BBSRC (to C. A. R.),the National Institutes of Health (AT006526 and HL077256 to Z. U.; AG022873 and AG025063 to S. N. A.), and the DFG Cluster of Excellence ‘‘Future Ocean’’ (to E. P.)

Outcomes associated with familial versus nonfamilial atrial fibrillation:a matched nationwide cohort study

BACKGROUND: We examined all‐cause mortality and long‐term thromboembolic risk (ischemic stroke, transient ischemic attack, systemic thromboembolism) in patients with and without familial atrial fibrillation (AF). METHODS AND RESULTS: Using Danish nationwide registry data, we identified all patients diagnosed with AF (1995–2012) and divided them into those with familial AF (having a first‐degree family member with a prior AF admission) and those with nonfamilial AF. We paired those with and without familial AF according to age, year of AF diagnosis, and sex in a 1:1 match. Using cumulative incidence and multivariable Cox models, we examined the risk of long‐term outcomes. We identified 8658 AF patients (4329 matched pairs) with and without familial AF. The median age was 50 years (interquartile range 43–54 years), and 21.4% were women. Compared with nonfamilial AF patients, those with familial AF had slightly less comorbid illness but similar overall CHA (2) DS (2)‐VASc score (P=0.155). Median follow‐up was 3.4 years (interquartile range 1.5–6.5 years). Patients with familial AF had risk of death and thromboembolism similar to those with nonfamilial AF (adjusted hazard ratio 0.91 [95% CI 0.79–1.04] for death and 0.90 [95% CI 0.71–1.14] for thromboembolism). CONCLUSIONS: Although family history of AF is associated with increased likelihood for development of AF, once AF developed, long‐term risks of death and thromboembolic complications were similar in familial and nonfamilial AF patients

Genetic variation and gene expression across multiple tissues and developmental stages in a non-human primate

By analyzing multitissue gene expression and genome-wide genetic variation data in samples from a vervet monkey pedigree, we generated a transcriptome resource and produced the first catalog of expression quantitative trait loci (eQTLs) in a nonhuman primate model. This catalog contains more genome-wide significant eQTLs per sample than comparable human resources and identifies sex- and age-related expression patterns. Findings include a master regulatory locus that likely has a role in immune function and a locus regulating hippocampal long noncoding RNAs (lncRNAs), whose expression correlates with hippocampal volume. This resource will facilitate genetic investigation of quantitative traits, including brain and behavioral phenotypes relevant to neuropsychiatric disorders

The Lick AGN Monitoring Project 2011: Dynamical Modeling of the Broad Line Region in Mrk 50

We present dynamical modeling of the broad line region (BLR) in the Seyfert 1 galaxy Mrk 50 using reverberation mapping data taken as part of the Lick AGN Monitoring Project (LAMP) 2011. We model the reverberation mapping data directly, constraining the geometry and kinematics of the BLR, as well as deriving a black hole mass estimate that does not depend on a normalizing factor or virial coefficient. We find that the geometry of the BLR in Mrk 50 is a nearly face-on thick disk, with a mean radius of 9.6(+1.2,-0.9) light days, a width of the BLR of 6.9(+1.2,-1.1) light days, and a disk opening angle of 25\pm10 degrees above the plane. We also constrain the inclination angle to be 9(+7,-5) degrees, close to face-on. Finally, the black hole mass of Mrk 50 is inferred to be log10(M(BH)/Msun) = 7.57(+0.44,-0.27). By comparison to the virial black hole mass estimate from traditional reverberation mapping analysis, we find the normalizing constant (virial coefficient) to be log10(f) = 0.78(+0.44,-0.27), consistent with the commonly adopted mean value of 0.74 based on aligning the M(BH)-{\sigma}* relation for AGN and quiescent galaxies. While our dynamical model includes the possibility of a net inflow or outflow in the BLR, we cannot distinguish between these two scenarios.Comment: Accepted for publication in ApJ. 8 pages, 6 figure

Effectiveness of moving on: an Australian designed generic self-management program for people with a chronic illness

Background: This paper presents the evaluation of “Moving On”, a generic self-management program for people with a chronic illness developed by Arthritis NSW. The program aims to help participants identify their need for behavior change and acquire the knowledge and skills to implement changes that promote their health and quality of life. Method: A prospective pragmatic randomised controlled trial involving two group programs in community settings: the intervention program (Moving On) and a control program (light physical activity). Participants were recruited by primary health care providers across the north-west region of metropolitan Sydney, Australia between June 2009 and October 2010. Patient outcomes were self-reported via pre- and post-program surveys completed at the time of enrolment and sixteen weeks after program commencement. Primary outcomes were change in self-efficacy (Self-efficacy for Managing Chronic Disease 6-Item Scale), self-management knowledge and behaviour and perceived health status (Self-Rated Health Scale and the Health Distress Scale). Results: A total of 388 patient referrals were received, of whom 250 (64.4%) enrolled in the study. Three patients withdrew prior to allocation. 25 block randomisations were performed by a statistician external to the research team: 123 patients were allocated to the intervention program and 124 were allocated to the control program. 97 (78.9%) of the intervention participants commenced their program. The overall attrition rate of 40.5% included withdrawals from the study and both programs. 24.4% of participants withdrew from the intervention program but not the study and 22.6% withdrew from the control program but not the study. A total of 62 patients completed the intervention program and follow-up evaluation survey and 77 patients completed the control program and follow- up evaluation survey. At 16 weeks follow-up there was no significant difference between intervention and control groups in self-efficacy; however, there was an increase in self-efficacy from baseline to follow-up for the intervention participants (t=−1.948, p=0.028). There were no significant differences in self-rated health or health distress scores between groups at follow-up, with both groups reporting a significant decrease in health distress scores. There was no significant difference between or within groups in self-management knowledge and stage of change of behaviours at follow-up. Intervention group attenders had significantly higher physical activity (t=−4.053, p=0.000) and nutrition scores (t=2.315, p= 0.01) at follow-up; however, these did not remain significant after adjustment for covariates. At follow-up, significantly more participants in the control group (20.8%) indicated that they did not have a self-management plan compared to those in the intervention group (8.8%) (X2=4.671, p=0.031). There were no significant changes in other self-management knowledge areas and behaviours after adjusting for covariates at follow-up. Conclusions: The study produced mixed findings. Differences between groups as allocated were diluted by the high proportion of patients not completing the program. Further monitoring and evaluation are needed of the impact and cost effectiveness of the program. Trial registration: Australian New Zealand Clinical Trials Registry: ACTRN1260900029821

Is childhood trauma associated with lifetime suicide attempts in women with bulimia nervosa?

The purpose of this study was to explore the association between specific forms of childhood abuse and neglect with lifetime suicide attempts in women with bulimia nervosa (BN)

Cultivation of common bacterial species and strains from human skin, oral, and gut microbiota.

BACKGROUND: Genomics-driven discoveries of microbial species have provided extraordinary insights into the biodiversity of human microbiota. In addition, a significant portion of genetic variation between microbiota exists at the subspecies, or strain, level. High-resolution genomics to investigate species- and strain-level diversity and mechanistic studies, however, rely on the availability of individual microbes from a complex microbial consortia. High-throughput approaches are needed to acquire and identify the significant species- and strain-level diversity present in the oral, skin, and gut microbiome. Here, we describe and validate a streamlined workflow for cultivating dominant bacterial species and strains from the skin, oral, and gut microbiota, informed by metagenomic sequencing, mass spectrometry, and strain profiling. RESULTS: Of total genera discovered by either metagenomic sequencing or culturomics, our cultivation pipeline recovered between 18.1-44.4% of total genera identified. These represented a high proportion of the community composition reconstructed with metagenomic sequencing, ranging from 66.2-95.8% of the relative abundance of the overall community. Fourier-Transform Infrared spectroscopy (FT-IR) was effective in differentiating genetically distinct strains compared with whole-genome sequencing, but was less effective as a proxy for genetic distance. CONCLUSIONS: Use of a streamlined set of conditions selected for cultivation of skin, oral, and gut microbiota facilitates recovery of dominant microbes and their strain variants from a relatively large sample set. FT-IR spectroscopy allows rapid differentiation of strain variants, but these differences are limited in recapitulating genetic distance. Our data highlights the strength of our cultivation and characterization pipeline, which is in throughput, comparisons with high-resolution genomic data, and rapid identification of strain variation